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We discuss how the baryon acoustic oscillation (BAO) signatures in the galaxy power spectrum 
can distinguish between modified gravity and the cosmological constant as the source of cosmic 
acceleration. To this end we consider a model characterized by a parameter n, which corresponds 
to the Dvali-Gabadadze-Porrati (DGP) model if n = 2 and reduces to the standard spatially flat 
qq 1 cosmological constant concordance model for n equal to infinity. We find that the different expan- 

' sion histories of the modified gravity models systematically shifts the peak positions of BAO. A 

, preliminary analysis using the current SDSS LRG sample indicates that the original DGP model 

■ is disfavored unless the matter density parameter exceeds 0.3. The constraints will be strongly 

tightened with future spectroscopic samples of galaxies at high redshifts. We demonstrate that 
\Q , WFMOS, in collaboration with other surveys such as Planck, will powerfully constrain modified 

f*"^ . gravity alternatives to dark energy as the explanation of cosmic acceleration. 
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£j ! I. INTRODUCTION 

Exploring the origin of the cosmic acceleration is one of the most challenging problems in cosmology. The nature 
of the cosmic expansion history can be investigated using a variety of cosmological observations: the anisotropics in 
the cosmic microwave background at z = 1100, the Hubble diagram of supernovae, cosmic shear statistics of galaxies 
with photometric redshifts, and spatial clustering of galaxies with spectroscopically measured redshifts. In the present 
paper, we focus on the fourth method which relies on extracting the baryon acoustic oscillations (BAO) in the high-z 
galaxy power spectrum Q, 0, 0, 0, IE E| ■ 

Recently the clear detection of the baryon oscillations has been reported @,EHUE3 already yielding some constraints 
surely on dark energy and particularly cosmic curvature. Based on this success, and the fact that the systematic errors 
in the BAO method appear easier to control (due to the fact that the characteristic scale of 100 Mpc is fundamentally 
linear), many BAO projects for the future are being discussed. For example, the AAOmega spectrograph on the 2dF 
system will be able to observe a large number of galaxies in the range of redshift 0.5 < z < 0.9 hll. The Fibre 
Multi Object Spectrograph (FMOS) on Subaru will be able to measure redshifts of more than 6 x 10 5 Ly-a emitting 
galaxies around z = 1.5 in near future 01 • The Large sky Area Multi- Object fiber Spectroscopic Telescope (LAMOST) 
project in China has a capability of performing large redshift-survey [13, while HETDEX (Hobby-Eberly Telescope 
Dark Energy Experiment) aims to measure redshifts of more than 10 6 Ly-a emitting galaxies using VIRUS (Visible 
Integral-field Replicable Unit Spectrograph) 13 • Furthermore, on a longer timescale, very large galaxy surveys might 
be performed with a space-based telescope [l5j and the SKA radio telescope |16j. 

In this paper we concentrate on the constraints that can be expected from the planned WFMOS. It is a BAO survey 
that aims to measure the redshift of more than two millions of high-z galaxies with a large field of view and multi-fiber 
spectrograph ( 4000 fibers) on a large ground-based telescope such as Subaru 0|. While the optimal geometry of 
WFMOS is not finalized yet we will use the fiducial baseline concept. One aim of the present paper is to demonstrate 
the power of WFMOS to extract the baryon oscillation in the galaxy power spectrum and the resulting tests that can 
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be made of modified gravity models which are putative explanations for the acceleration of the cosmos. 

Although the cosmological constant is the simplest model leading to acceleration, new dynamical degrees of freedom 
(e.g., a scalar field) with effective negative pressure may be the source of acceleration. An interesting possibility is 
that acceleration results not from an additional repulsive force but through the weakening of traditional Einstein 
gravity on very large scales. The archetypal model in this class is the Dvali-Gabadadze-Porrati (DGP) model which 
was developed in the context of the brane world scenario, which includes a mechanism to explain self-acceleration 
in the late time universe [l8l lisj . The DGP model and variants [20L I2H I22I I23I ] will serve as a testbed for testing 
the ability of WFMOS to distinguish them from a cosmological constant. We note here that the DGP model is not 
without subtleties and problems at the quantum level as discussed in 0, HE H3 ■ 

Very recently, the evolution of the large-scale cosm olog ical perturbations of the DGP cosmological model have 
been studied in detail with rather general assumptions [23 compared to earlier works |H, |2j| |3(J • Koyama has also 
presented a generalized model that interpolates between the DGP model and the CDM model with a cosmological 
constant (ACDM) in general relativity [31|. We will see that WFMOS may well provide not only a chance to test 
dark energy models but also to put general relativity to severe tests on large scales. To achieve this we investigate 
the theoretical predictions of modified gravity models in the galaxy power spectrum from a survey with WFMOS. 

So far, there are a few works which focus on the nature of the large scale structure in a modified gravity model 
|32l 133. l34j. In the present paper, we consider the DGP cosmological model and its variants from the point of view 
of testing the gravity theory with the baryon oscillation, and investigate a future prospect assuming the WFMOS 
project for the first time. We also investigate the current consistency with observation utilizing the power spectrum. 
A similar problem is considered in [U US 13 , however, these works make use of a result on the distance from the 
SDSS LRG correlation function analysis in 

This paper is organized as follows: In section 2, a brief review of the DGP cosmological model and its variants is 
presented. In section 3, our theoretical modeling of the power spectrum and our basic statistic is explained. Then, 
a current comparison of the theoretical model with the observational result with the luminous red galaxy (LRG) 
sample in the SDSS is presented in section 4. In section 5, we discuss the baryon oscillation features, assuming a 
future WFMOS sample, as a test of the modified gravity model. We discusses how one can improve the constraints 
on modified gravity models from future WFMOS samples. The last section is devoted to a summary and conclusions. 
Throughout this paper, we use units in which the speed of light and the Planck constant are unity, c = h = 1. We 
adopt the Hubble parameter: Hq = 72 km/s/Mpc, unless we mention explicitly. 



II. THE DGP COSMOLOGICAL MODEL AND VARIANTS 

We start with a brief review of the DGP cosmological model (see e.g., H3)- The DGP model is based on the brane 
world scenario, and is constructed by embedding a (3+l)dimensional brane in a (4+1) dimensional bulk with infinite 
volume, with action 

S = ~ I d 5 X^^)R^ - ^ / d^x^g^R^ + / <£x^^)L m + 5 g .h., (2.1) 



16-7T J ' 167T 

where M, (M p \) is the fundamental Planck mass in the 5-dim (4-dim) space-time, g^ (g^) and (R^) denote 
the determinant and the Ricci scalar of the 5-dim (4-dim) metric, respectively, and C m is the matter Lagrangian 
on the brane. The final term, the Hawking-Gibbons term Sg.h. 1 is added so as to reproduce the appropriate field 
equation in a space-time with a boundary. The term with i?' 4 ' in the above action is assumed to be induced by 
quantum effects in the matter sector on the brane. 

The crucial cross-over scale, r c , is defined by the ratio of the Planck scales 

This characterizes the region where gravity switches from being 4-dimensional (scales less than r c ) to being 5- 
dimcnsional (scales greater than r c ) which causes a modification of the laws of gravity on cosmologically larger 
scales if r c ~ H^ 1 . As a result so-called self-acceleration may appear. The cosmological solution of the DGP model 
give 



2 H _ 8tt P 



where Hit) is the Hubble expansion rate, pit) is the matter density, and M pl 2 is regarded as the 4-dimensional 
gravitational constant, G. 



3 



Cosmological perturbations have also been studied in the modified DGP-like models |3jJ, originally proposed by 
Dvali and Turner j3]j ■ These variants are phenomenological models interpolating between the DGP model and ACDM. 
Adopting this modified DGP-like model, the modified Friedmann equation is 

where r c and n are the parameters of the model. In this paper, we assume spatially flat sections, implying that the 
parameter r c is related to the cosmological parameters by (Hor c ) 2 ' n ~ 2 = 1 — fl m , where fl m is the matter density 
parameter. 

This model is the DGP model for n = 2, and it reduces to ACDM for n — > oo. Following the analysis of the modified 
DGP-like model [3lj. the evolution equation for the linear density perturbation (linear growth factor) obeys 

D 1+ 2HD 1 = ^(l + ±y Dl , (2.5) 

where 

(3 = 1- n{Hr c ) 2{n -^/ n (l + 2(n - l)H/3nH 2 ^j . (2.6) 

In the present paper, we test these modified DGP models by extracting the baryon oscillations of simulated and real 
galaxy power spectra. 

Concerning the linear growth factor , th e time evolution slightly depends on the treatment and assumptions of the 
perturbation equation (23, |28|, H^, LM H3 . However, the difference is small, and the effect on the baryon oscillation 
feature is negligible in the cases studied in the present paper. In a strict sense, the growth factor is the expression 
under the sub-horizon approximation. Hence the growth factor may depend on wavenumber at large scales, which 
could lead to additional scale dependence [38l| . 



III. BARYON ACOUSTIC OSCILLATION IN THE MODIFIED DGP MODEL: THEORETICAL 

PREDICTIONS IN LINEAR THEORY 



The reason why clustering statistics are sensitive to the expansion history of the universe comes primarily from the 
baryon acoustic oscillations (BAO), which imprints a characteristic scale on the galaxy distribution that acts as a 
standard ruler, which we briefly summarize here (see also e.g. US USUI!). The origin of the BAO in the matter 
power spectrum can be understood as the velocity fluctuation of the baryon fluid at the decoupling time [3^, |4jj . 
The characteristic scale of the baryon oscillation is determined by the sound horizon at decoupling, which depends 
on the total matter density and baryon densities. Because we can measure this scale in both the transverse and 
radial directions the BAO yields both the angular diameter distance and Hubble parameter at that redshift. Thus, 
the precise measurement of the BAO scale from galaxy power spectrum can put important constraints on the cosmic 
expansion history. 

Our model predictions of galaxy power spectra in the modified DGP model proceed as follows (see |4lL |42| for 
details). In the linear regime, the real-space mass power spectrum P m &ss(q', z) is translated to the galaxy power 
spectrum in redshift-space as 



■fgal(?lh?x;«) = b 2 {z) 



1 din D\(z) ( gu 
b(z) dlna(z) \ q 



Pmass \Q'i Z) 



(3.1) 



where a{z) is the scale factor at z, and q ± is the parallel and perpendicular component of the comoving wavenumber 
to the line-of-sight direction. We assume a linear scale-independent galaxy bias, b(z). 

Mapping the length scale from the observed distribution, i.e.. redshift and angular separations Sz and 68, becomes 
sensitive to the cosmological parameters at high redshift s[43l El EEj . If one adopts a fiducial set of cosmological 
parameters in estimating the power spectrum from the observed galaxy distribution, the result is necessarily distorted 
from the actual spectrum, but the effect can be modeled theoretically. Suppose that the Hubble parameter in the 
fiducial model is H^\z) instead of the true one, H(z), which can be obtained by numerically solving the modified 
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Friedman equation l|2.4|) . Then, comoving radial distances in the fiducial and true universes are given by 1 



r (f) (z) 



dz' 



r{z) 



dz' 



o H(z'Y 



/o H(t)(z<y 

The parallel and perpendicular comoving distances that correspond to the observables 5z and 59, respectively, are 

1 



(3.2) 



dr(z) 

Xii = ; OZ 



-8z, x x = r(z)S9. 



(3.3) 



dz H(z) 

Therefore their counterparts in the model universe with the adopted fiducial parameters are related to the above as 



(f) _ H(z) 
X " HV)(z) 

The corresponding relations for the wavenumbers read 



(f) _ ifHfL 

r(z) 



which yield 



= g|l 2 + ^ 2 = [fc( f )] 2 



/r 



"( f )( z ) 



+ (1-M 2 ) 



r (f)( z ) 



Hit){z 

Then, the power spectrum in the fiducial universe is computed from 

P^(k^,fl,z)d 2 k { ! ) dklP = P gal (q„q ± ;z)d 2 q ± dq 
Combining equations (|3.1|) . I|3.5|) . and (|3.(jl) . one obtains 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



F (f) (fc (f) ,Ai,z) = b 2 (z) 



1 dlnL>!(z) 
6(z) dlna(^) 



M 2 + (l-M 2 ) 



fH^(z)r^(z) 



H(z)r(z) 



H(z) 
H( { \z) 



-(f) 



-(z) 



X P„ 



ff(f) 

i?( f )(z) 



+ (1-M 2 ) 



■( f )(z; 

r(z) 



(3.8) 



Finally, the estimated power spectrum is given by integrating equation 13.81) over the light-cone: 



P(k) 



d/j, 



dz dr(i ^ Z K ^(z) 2 fi(z) 2 i;(z, k) 2 P^(k, fi] z) 



dz 



dr(t){z) r^(z) 2 n(z) 2 ^(z,k) 2 



(3.9) 



dz 



where we write k^ simply as k. The mean number density of galaxies in the cosmological redshift space, n(z), is not 
directly observable, but rather inferred from the observed number count dN(z)/dz: 



dz ~ 1 ' HW(z) ' 



(3.10) 



We may introduce a weight factor, ip(z, k), which is set as 1/[1 + P(k, /i, z)n(z)] in the optimal weighting scheme. In 
the present paper, we adopt ip(z,k) — l/n(z) for the SDSS LRG sample following ref. 9], and tp(z,k) — 1 for future 
WFMOS samples for simplicity. 



1 The comoving distance in the fiducial universe is denoted by s(z) in the previous paper 113 . instead of A l \z) 
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The analysis presented in this paper focuses on the local power-law index of the power spectrum: 

d\nP(k) _ 1 dP(k) 
d\nk ~ P{k) d\nk ' 



(3.11) 



In general, d In P(k)/d In k does not contain all the information encoded in P(k) so constraints from d\nP(k)/d\nk 
will be somewhat weaker than those from P(k). However, d\nP{k)/dh\k may be useful in extracting the baryon 
oscillation feature because, as long as the oscillations are clearly detected the principal part of the constraints on 
the cosmic expansion history using the power spectrum come from the oscillations |4q . In addition, d\nP(k)/d\nk 
does not depend on the amplitude of the power spectrum. Thus, it is not sensitive to the bias and the growth factor 
which provides a desirable level of model-independence. In our theoretical predictions, we obtain d In P(k)/d In k by 
differentiating P(k) with respect to k which does introduce some level of additional noise. 

The fiducial cosmological model that we choose in the following analysis is a spatially-flat A CDM with Q$ — 0.27, 
fiW = 1 - 0$, nb _ 0.044, and h = 0.72. Therefore 



ff (f) 0) = H Q ^n$(l + z)3 + l-ni\ (3.12) 



Since we assume scale- independent linear growth factor in the modified DGP model (ea. |2.5j ). P m ass(<?; z), except 
for its overall amplitude, is exactly the same as that in the fiducial ACDM model. The spectrum is computed u sing 
the linear transfer function by Eisenstein and Hu [3^, and the primordial spectral index n s = 0.95 and o% = 0.8 |47|. 
For definiteness we adopt ip(k, z) = 1, and the Fry's bias model: 

b{ z ) = l + b ±zl (3.13) 
D 1 {z) 

with b = 1.5 for the samples WFMOS1 and WFMOS2, and with b = 2 for the sample SDSS LRG (see below for 
details). The reasonable change in the bias model does not alter our results for d\nP(k)/d\nk. However the error 
does depend on the bias amplitude. 

Figs. 1 and 2 plot r(z)/r^\z) and H^\z)/H(z), respectively, as function z for the original DGP model (n = 2; 
the dashed curve), and the modified DGP models with n = 4 (dotted curve) and n = 8 (dash-dotted curve). Here 
the density parameter is fixed as fl m = 0.27. The solid curve is the ACDM model, which is trivially r{z)/r^\z) = 
H^\z)/H(z) = 1, because it is equivalent to the fiducial model. For clearly differentiating models, the deviation 
of r(z)/r(^(z) and H^\z)/H(z) from unity is important. Therefore, the sample of the redshift larger than 0.5 is 
effective for selecting the models. 

In the rest of this section, we consider theoretical predictions of the modified gravity model for d In P(k)/d In k. 
Here we assume a future fiducial WFMOS survey in the range of redshifts 0.5 < z < 1.3 (z m i n = 0.5 and z max = 1.3; 
see section 5 for details) . Fig. 3 shows the theoretical curves of d In Pjd In k for the original DGP model (n — 2, dashed 
curve), the DGP-like models with n = 4 and n — 8 (dotted and dashed-dotted curves) and for the ACDM model (solid 
curve, the fiducial model). We have fixed the spectral index of the initial spectrum to the value n s = 0.95 and the 
density parameters to fib = 0.044, fl m = 0.27. The difference of predictions between the fiducial A CDM model and 
the modified gravity models entirely comes from the geometric distortion, i.e., the model-dependent relation between 
q and k, equation (|3.5|) . The apparent geometric distortion disappears only when the actual universe has the same 
H^(z) and r^\z) for z < ^max of the fiducial model. 

It is instructive to consider how the theoretical predictions are sensitive to the cosmological parameters in order to 
understand the degeneracy among different parameters. For this purpose, we plot the results for different A CDM 
models by varying the values of fi m and fib in Figs. 4 and 5, while the other parameters are the same as the fiducial 
A CDM model in Fig. 3. Note again that the wavenumber k in Figs. 4 and 5 is also computed under the assumption 
of the fiducial model parameters and should be regarded as k^\ 

Figs. 4 and 5 indicate that varying fi m and fib in spatially-flat CDM models mainly affects the amplitude of the 
BAO, but the induced shift of the peak positions is weak, at least much weaker than the difference between the DGP 
model and the fiducial model illustrated in Fig. 3. This also applies to the change of the spectral index, n s . Therefore 
we conclude that the degeneracy among those parameters space is not so crucial, although it is clear that a prior 
constraint on these parameters is important and complementary to the BAO test. 



IV. CURRENT CONSTRAINTS: SDSS LRG SAMPLE 



In this section, we utilize the theoretical framework developed in the previous section to place constraints on modified 
DGP models from existing large-scale structure data. For this purpose, we use the published power spectrum of galaxy 
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clustering by Hiitsi Q based on 51763 LRGs from Data Release Four (DR4) of the Sloan Digital Sky Survey (SDSS). 
In this work, they present a discovery of the BAO in the power spectrum and compute the covariance matrix as well, 
which is ideal to apply our current methodology. 

We fit our cosmological models to the power spectrum data by varying Q m and assuming spatial flatness both 
for the ACDM and modified DGP models. In Figure 6, we present the best fit A CDM model (solid line) and the best 
fit DGP model (dashed line) to the SDSS LRG power spectrum of [9j. For comparison, we also present the predicted 
DGP model (dotted line) with the same cosmological parameters as the best fit A CDM model, i.e., with fi TO = 0.32 
and fib = 0.045. For these fits, we have used z min = 0.16 and z max = 0.47, as well as determined the comoving 
number density n(z) (see Fig. 4 of 0) in the fiducial model cosmology (assuming the first year WMAP parameters). 
The error bars were computed via Monte Carlo simulations and the covariance matrix provided by 9]. To minimize 
non-linear effects, we only fit to the SDSS LRG power spectrum at scales of k < 0.2 ftMpc -1 . When quasi-nonlinear 
effects properly taken accounted for[^[4£|,[4£j,[5(j, we may increase the range of scales used in this analysis and thus 
improve the statistics. 

In Figure 7 (a), we show the contours of A\ 2 in the fl m versus Qb/ttm parameter plane with respect to the best fit 
ACDM and DGP models shown in Figure 6. We compute the x 2 using 

\d\nP/d\nk\° h \ -dlnP/dlnfc|* he °l 2 
Y [A(d\nP/dlnk)} 2 ' K ' ' 

where d\nP/d\nk\° hs and A(d\nP/d\nk) are the measured value and the errors in Figure 6 respectively. The 
way of finding A(d\nP/d\nk) is similar to that in the next section. We generate mock P(k) data from Monte 
Carlo simulations with the covariance matrix in Q . Then, we evaluated the variance of d In P/d In k (see also next 

section). Also, din P/<i In /c| tho ° is the corresponding theoretical value and is evaluated from P(k) thco sampled at the 
wavenumber bin interval of Ak = 0.02 /iMpc" 1 as in the observational data 0. The other cosmological parameters 
are held at the values set by our fiducial model. 

In Figure 7(a), we plot the contour levels of A\ 2 = 2.3 (inner curve) and 6.2 (outer curve), which correspond to the 
one sigma and two sigma confidence levels for the x 2 distribution (with two degrees-of- freedom). Clearly, the DGP 
model favors a higher value of Q m than the ACDM model (see ref. J23|). Figure 7 (b) is the same as Fig. 7 (a) but 
with assuming Hq — 66 km/s/Mpc. Thus, the alternation of the Hubble parameter does not change our conclusion. 

In figure 8, we show the contours of A\ 2 for the modified DGP models in the fl m versus 1/n parameter plane. For 
reference, n = 2 corresponds to the original DGP model, while n = oo is the ACDM model. 

Figures 7 and 8 show the potential for differentiating between DGP models and the standard ACDM model if 
independent measurements of tt m , or the other cosmological parameters, can be obtained from additional data. We 
do not attempt such an analysis here because of the small redshift range of the SDSS LRG sample that results in a 
relatively small predicted difference in the location of the BAO peaks between the DGP models and ACDM model. 
Furthermore, we have not included non-linear effects. Instead, we look forward to the next generation of galaxy 
redshift surveys where the redshift baseline will be much larger and therefore, the expected differences between the 
DGP and ACDM models will be greater. 

Limiting to the ACDM model, we briefly compare our result with that in other analyses. In Hiitsi investigated 
cosmological constraints from the SDSS LRG power spectrum. Our result on fl m is larger than his result. This 
might come from the difference in the treatment of the non-linear effect. In addition, his result is obtained by being 
combined with the WMAP result, then the comparison is not straightforward. Our result can be rather consistent 
with that in 1521. 



V. FUTURE CONSTRAINTS: WFMOS SAMPLES 



As stated above, we now consider the expected constraints on DGP models from future galaxy samples at higher 
redshifts. In particular, we consider the two galaxy redshift samples proposed by the WFMOS experiment; "WF- 
MOS1", which will contain 2.1 x 10 6 galaxies, over 2000 deg 2 , at 0.5 < z < 1.3 and "WFMOS2" , which will contain 
5.5 x 10 5 galaxies, over 300 deg 2 , at 2.3 < z < 3.3. The corresponding comoving mean number density for these two 
samples is n = 5 x 10~ 4 (h~ 1 Mpc)~ 3 and n = 4x 10~ 4 (h~ 1 Mpc)~ 3 , respectively. 

We determine the likely error on the power spectrum of these two WFMOS galaxy samples using Monte Carlo 
simulations. In detail, the error on P(k) is obtained by the analytic formula of Feldman ct al. 53] . More specifically, 
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we use (see also |42],|5J) 



dp I dz 1 ' r^(z) 2 n(z) 4 iP{k,z)' l [P(k,fi;z) + l/n} 2 

= 2^- - ^ . ' , m „ , (5.D 



Aft 



d«^^r< f >(«) 2 n(z)V(M) 2 ' 



where AV4(= 47rfc 2 Afc) denotes the volume of the shell in the Fourier space and Aft is the survey area. Note that 
we adopt ip(k,z) = 1 here. We also create mock data from Monte Carlo simulations. We first generate P(ki) with 
random errors from a Gaussian distribution, where fej = iAk. Then, we evaluated <i In P/<i In fc at the wavenumber 
(hi + k i+ i)/2 from the nearest two bins by (P(k i+ i) — P(fcj))(fcj + k i+ i)/(P(ki) + P(k i+ i))/Ak. The quoted error 
bars are one sigma derived from 10 5 mock realizations. 

In Figure 9, we present theoretical predictions of d\nP/d\nk for the WFMOS1 galaxy sample based on the ACDM 
(solid line) and modified DGP model (dashed line) presented in Figure 3. The data points with error bars are 
derived from our mock data with the wavenumber binsizc of Ak = O.OlfcMpc -1 . We then compare the predictions 
of the ACDM model and the DGP model (n — 2) with the mock data in Figure 9 using data points in the range 
0.02 < k < 0.2 /iMpc -1 . Based on chi-square, we find that these two models are different by 6 a. The statistical 
difference, with respect to the ACDM model, decreases to 2a for a modified DGP model with n = 3, while modified 
gravity models with n > 4 are only statistically different at the la level. Therefore, Figure 9 demonstrates that the 
BAO analysis of the WFMOS1 galaxy sample will be able to discriminate between the standard ACDM cosmology 
and the original DGP models. Moreover, such analyses one will put stringent constraints on the modified DGP models 
with n > 3 especially when combined other independent cosmological observations. 

For the actual analysis of the WFMOS galaxy power spectra, the binsize (Afc) will be determined to ensure the 
statistical independence of adjacent bins, which will depend on the final survey geometry and the number density of 
observed galaxies. To test the effect of such issues on our results, we present in Figure 10 new predictions based on a 
binsize of Afc = 0.02/iMpc . Furthermore, in Figure 11, we adopt the nearest three points in computing dP/dk from 
the mock data. Taken together, Figures 10 and 11 demonstrate our results are robust against the details of binning 
the data. Nevertheless the statistical significance of the constraints should be carefully examined by properly taking 
account of the covariance matrix analysis of different data points, which is beyond the scope of the present paper. 

Figure 12 is similar to Figure 9 but now for the higher redshift WFMOS2 galaxy sample discussed above. The 
theoretical curves are almost the same as those for the WFMOS 1 sample, but the error bars are now a factor 2 to 3 
larger, because the theoretical curves are determined by the factors r(z) /r^\z) and H^\z) / H(z) (Figures 1 and 2) 
which are almost constant over the redshift range 0.5 < z < 3. 

To estimate future cosmological constraints, Figure 13 shows the contour of A\ 2 in the fl m and 1/n parameter plane 
for both the WFMOS1 (solid curve) and WFMOS2 (dashed curve) galaxy samples. This is similar to Figure 8 for the 
existing SDSS LRG sample. In this figure, the underlying cosmological model is ACDM with Sl m = 0.27, n s = 0.95 and 
fib = 0.044. Using Monte Carlo simulations, we create 100 mock datasets and computed d\nP(k)/d\nk with a binsizc 
of Afc = 0.01/iMpc -1 (Figure 9) and then we compute the average A% 2 statistic. The errors on the power spectrum 
were determined from the effective volume of the WFMOS1 and WFMOS2 samples, and the amplitude of the power 
spectrum (the volume of the WFMOS2 sample is assumed to be a quarter of WFMOS1 volume). Furthermore, we 
assumed the same bias parameter (bo = 1.5) which causes a decrease in the amplitude of power spectrum of WFMOS2 
by a factor 2. 

We also consider here a new WFMOS sample of galaxies (called "WFMOS3" ) which consists of 2.2 x I0 6 galaxies 
over 1200 deg 2 at 2.3 < z < 3.3. This gives the same number density as the WFMOS2 sample, but with almost same 
volume as the WFMOS1 sample. We also assume bo = 1.9, so that the amplitude of the power spectrum is almost 
same as the WFMOS1 sample. We present in Figure 13 the contour of A% 2 for the WFMOS3 sample (dotted curve), 
which demonstrates that such a sample would provide similar cosmological constraints as the WFMOS 1 sample. 

The error of d In P/d In fc can also be improved by increasing the mean number density of galaxies. The range of the 
relevant wavenumber (the maximum wavenumber) is also important for constraining the parameters. These factors 
will be able to improve the constraint. 

Finally, we stress that our analysis is based on the fitting of din P(fc)/<i In fc. In general, however, the fitting of P(k) 
gives more stringent constraint. Then, the constraint in this paper can be improved when using the fitting of P(k) 
directly. But, in our analysis, we fixed the cosmological parameters n s and Uncertainty of these factors might 
weaken the constraint. 
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VI. DISCUSSION AND CONCLUSIONS 



We have considered here how models of modified gravity affect the BAO signature in the power spectra of galaxy 
clustering both for current and future galaxy samples. As a specific model of the modified gravity, we adopt the 
modified DGP model[3l| which is characterized by an index n; the model with n = 2 corresponds to the original 
DGP model^li while n = oo approaches the concordance, spatially flat, ACDM model in general relativity. We have 
shown that the different expansion history in the modified gravity models shifts the peak positions of oscillations 
relative to the ACDM model. These predicted shifts in the BAO can potentially be used to distinguish between the 
ACDM model and the modified DGP model. 

We tested our predictions for the BAO against existing data from the SDSS LRG sample and found we can already 
rule out the original DGP modified gravity model (with n = 2) if we assume fi m < 0.3 with the flat geometry, although 
we have not carried out a full analysis of the parameter space (cf. [28l Isfil]). For future WFMOS galaxy samples, 
there is great promise for putting stringent constraints on the modified gravity. 

We also note that our work is based on the linear theory of density perturbation. The non-linear nature of the 
gravity force in the DGP cosmology is not clear at present. It might be expected that the modification of the Poisson 
equation would vanish by the non-linear nature. However, this problem remains unsolved. 

In our investigation, we focused on the quantity dlnP(fc)/dlnfc in extracting the oscillation feature. As demon- 
strated in the last section, the theoretical curve and the result of the simulation is slightly different. This comes 
from the finite binsizc effect, which becomes problematic when the bin of k is large. This gap might be resolved by 
introducing a new algorithm for estimating in d P/d\ n k. Following the conventional method of evaluating the power 
spectrum with a discrete number density field |42T . l53tl54| . the power spectrum can be obtained by using the estimator, 
adopting the constant weight factor, 



P(k) = / dsn(s) 



N 

E 



N/c 



e lk Sl - a J2 e' k Sj 



shot noise term, 



(6.1) 



where s, (1 < i < N) is the i-th galaxy's position of a real catalog, Sj (1 < j < N/a) is the j-th galaxy's position of 
a random catalog, and a a parameter chosen as a < 1. In the similar way, we can evaluate the differentiation of the 
power spectrum by introducing the estimator 



dP{k) 

' dk 



dsn(s) 2 I ik • Sjj 



JY 



N/a 



E ik ' s ii 



N 

E 



N/a 



-complex conjugate 



(6.2) 



This expression gives an alternative way to evaluate dP(k)/d\nk. The effectiveness of the use of this estimator will 
be investigated as a future work. 

Finally, we note the degeneracy between the DGP models used in this paper and other dark energy models. As 
demonstrated in |55l l56|. the expansion history of the DGP model can be reproduced by a dark energy model with 
an effective equation of state of w c r(z) ~ wo + w a x z/(l + z). The behavior of the galaxy power spectrum of the 
DGP model presented in this paper can thus be re-interpreted using this effective equation of state. This is because 
the WFMOS surveys only constrain the expansion history of the Universe and therefore, the measurement of the 
power spectrum alone can not differentiate between the DGP model and more complex dark energy models if they 
have the same expansion history. Observationally, weak lensing measurements can be more sensitive to the growth 
factor than the BAO. Therefore HypcrSuprime and other large weak lensing surveys will be ver y u seful in breaking 
this degeneracy while weak lensing in turn suffers from an fix degeneracy that BAO can break 57]. Together BAO 
and weak lensing offer a powerful synergy in the hunt for the origin of cosmic acceleration. We hope to return to this 
point in future. 
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FIG. 1: r(z)/r^\z) as function of z for the DGP model (dashed red curve), the DGP-like modified gravity models n = 4 
(dotted blue curve), n = 8 (dash-dotted green curve), and the ACDM model (solid black curve). Here the density parameter 
is fixed as Q™ = 0.27. 
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FIG. 2: Same as Fig. 1 but for H if) (z)/H(z). The DGP model (dashed red curve), the DGP-like modified gravity models 
n — 4 (dotted blue curve), n = 8 (dash-dotted green curve), and the ACDM model (solid black curve). Here Q m = 0.27 as 
Fig. 1. 



[57] G. Bernstein Astrophys. J. 637, 598 (2006) 
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FIG. 3: Theoretical curves of din P/dlnfc for the DGP model (dashed red curve), the DGP-like modified gravity models n = 4 
(dotted blue curve), n = 8 (dash-dotted green curve), and the ACDM model (solid black curve) assuming the sample WFMOS1. 
The modified gravity with large n approaches to the ACDM model. Here we adopted the initial power spectral index, n s — 0.95, 
and Q b = 0.044 and f2 m = 0.27. 




FIG. 4: Theoretical curves of din P/d In k in the ACDM model showing the effect of changing fl m : the solid curve has £l m = 0.27, 
the dashed curve has f2 m = 0.22, the dotted curve has Q m = 0.32. The other parameters are the same as those of Fig. 3, 
n a = 0.95 and Q b = 0.044. 
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FIG. 5: Theoretical curves of d In P/d In k in the ACDM model showing the effect of varying Qb'- the solid curve has Qb = 0.044, 
the dashed curve has Qb = 0.04, and the dotted curve has Q m = 0.05. The other parameters are the same as those of Fig. 3, 
n s = 0.95 and Q m = 0.27. 
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FIG. 6: din P/d In k of the SDSS LRG sample and the corresponding theoretical curves of the DGP model (dashed red curves) 
and the A CDM model (solid black curve). The thick curves are the best fit model which has fl m — 0.32 and Qb = 0.045 for 
the A CDM model, and Q m = 0.35 and Qb — 0.053 for the DGP model. The thin dashed curve is the DGP model with the 
same parameter as the A CDM model (f2 m = 0.32 and fij, = 0.045). The other cosmological parameter n 3 — 0.95 is the same 
as in Fig. 3. The error bar is evaluated using the simple simulation with the measured power spectrum and the covariance 
matrix in ref. [9]. 
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FIG. 7: (a) Contour of Ax 2 on the parameter space Sl m and Slb/Sl m - The other parameters are fixed as the same as those in 
Fig. 3. The solid curve is the A CDM model and the dashed curve is the DGP model. The contour level is A\ 2 = 2.3 (la; 
inner curve) and A\ 2 = 6.2 (2er; outer curve). Here n s — 0.95 as before, (b) Same as (a) but with h — 0.66. 




FIG. 8: Contour of Ax 2 on the parameter space Q m and 1/n in the DGP-like model, n = 2 corresponds to the DGP model. 
The other parameters are fixed as n s = 0.95 and Sib = 0.044. The contour level is \ 2 = 2.3 (la; inner curve) and \ 2 = 6.2 (2a; 
outer curve). 
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FIG. 9: Theoretical predictions for dlnP/dlnk assuming the sample WFMOS1 (corresponding to z < 1.3). The squares with 
error bar are evaluated with the simple simulation of the power spectrum for the A CDM model. The asterisks are the DGP 
model, but the error bar, which is almost the same as that of the A CDM model, is omitted for simplicity. Theoretical curves 
are the DGP model (dashed red curve) and the A CDM model (solid black curve). The parameters are the same as those in 
Fig. 3, i.e, n 3 = 0.95, fij, = 0.044, Q m = 0.27 and bo — 1.5. Here we used the nearest two points in the differentiation and 
Ak = 0.01 hMpc -1 . 
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FIG. 11: Same as Fig. 9, but for the case of the differentiation with the nearest three points. The parameters are chosen as 
n s = 0.95, fl b = 0.044, fl m = 0.27 and b = 1.5. 
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FIG. 12: Same as Fig. 9, but for the sample WFMOS2, which assumes the range of the redshift 2.3 < z < 3.3 with n = 
4 x 10 _4 (ft _1 Mpc) -3 , and a survey area 300 deg. 2 , which yields a total number 5.5 x 10 5 galaxies. The parameters are chosen 
as n s = 0.95, Q b = 0.044, Q. m = 0.27 and 6 = 1.5. 
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FIG. 13: (a, solid curve) Contour of A\ 2 on the parameter space Q m and 1/n in the DGP-like model, assuming the sample 
WFMOS1, which consists of 2.1 x 10 6 galaxies over 2000 deg 2 at 0.5 < z < 1.3 with the mean number density n = 5 x 
10~ 4 (ft _1 Mpc) -3 . The target model is the A CDM model, i.e., 1/n = and Q m = 0.27, and assumed the bias parameter 
bo = 1.5. Here we fixed n a = 0.95, Q.b = 0.044. The contour level is \ 2 — 6.2 (2a). (b, dashed curve) Same as (a) but 
for the sample WFMOS2, which consists of 5.5 x 10 galaxies over 300 deg 2 at 2.3 < z < 3.3 with the mean number density 
n = 4x 10~ 4 (fe _1 Mpc) -3 , and assumed the bias parameter bo = 1.5. (c, dotted curve) Same as (a) but for the sample WFMOS3, 
which consists of the 2.2 x 10 6 galaxies over 1200 deg 2 at 2.3 < z < 3.3 with the mean number density n — 4x 10 -4 (/i -1 Mpc)~ 3 , 
and assumed the bias parameter bo = 1.9. 



